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We consider propagation of sound pulses along the long axis of a Bose-Einstein condensed cloud 
in a very anisotropic trap. In the linear regime, we demonstrate that the square of the velocity of 
propagation is given by the square of the local sound velocity, (? = nUo/m, averaged over the cross 
section of the cloud. We also carry out calculations in the nonlinear regime, and determine how the 
speed of the pulse depends on its amplitude. 
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In a beautiful recent experiment, Andrews et al. have 
reported measurements of the propagation of sound 
pulses along the axis of an elongated cloud of Bose- 
Einstein condensed Na atoms Q. At the center of the 
trap a disturbance was created by turning a laser on 
or off. Due to the energy shift caused by the interac- 
tion of atoms with the laser field, this led to a change 
in the potential experienced by atoms, and thereby to a 
density disturbance, whose time development was stud- 
ied by phase-contrast imaging. In a dilute, homogeneous 
Bose gas at zero temperature, the speed of sound, cq, is 
given by Cq = dP/dp = noUo/m, where P is the pres- 
sure, p = mn is the mass density, m is the atomic mass, 
n is the density of atoms, and Uq = ^Trh^a/m is the 
strength of the effective two-body interaction, with a be- 
ing the scattering length for atom-atom collisions 
The experimental data seemed to be in accord with this 
result, if the density was taken to be that at the center 
of the cloud. Within the accuracy of the measurements, 
no evidence for a dependence of propagation speed on 
amplitude was seen. 

The purpose of this Letter is to investigate how propa- 
gation of a sound pulse is influenced by the inhomogene- 
ity of the cloud, and by nonlinear effects. The important 
result of our calculations is that, in the linear regime, the 
sound speed is given by the expression for the homoge- 
neous gas, but evaluated for the average density over the 
cross section of the cloud perpendicular to the direction 
of propagation. A second result is that nonlinear effects 
are substantial, and can give rise to changes in the ve- 
locity of the pulse of order of 10%, for relative density 
perturbations of order 15%. 

Consider a Bose-Einstein condensed cloud whose trans- 
verse dimensions (~ Rs_) are much less than the axial 
dimension (^^ Z). We study the propagation along the 
axis of the trap, of pulses whose characteristic spatial 
scale, /, in the axial direction is large compared with 
but small compared with Z. That this condition 
is satisfied in the experiments may easily be seen by in- 
spection of the density profiles in Fig. 2 of Ref. For 
the moment we make no detailed assumptions about the 



confining potential perpendicular to the long direction of 
the trap, which we take to be the z-axis, apart from the 
fact that its variations with respect to z should be small 
on a length scale I. The fact that the transverse dimen- 
sions of the cloud are small compared with the scale of 
the pulse, means that the density profile across the trap 
may be assumed to have its equilibrium form appropri- 
ate to the local number of particles per unit length. This 
follows from the fact that the timescale for adjustment 
of the profile ('^ R±/c), where c is a typical sound speed, 
is short compared with the time for passage of the pulse, 
l/c. The problem then becomes a one-dimensional one, 
and the sound pulse may be specified in terms of a lo- 
cal velocity, v{z), and a local density of particles, cr(z), 
per unit length. The latter is given in terms of the par- 
ticle density, n{x,y,z), where x and y are coordinates 
perpendicular to the axis of the trap, by 



a{z) ~ / dxdy n{x , y , z) . 
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The equations governing the motion of the condensate 
may be derived from the Gross-Pitaevskii equation, and 
they are just the equations of perfect fluid hydrodynam- 
ics, with the addition of a "quantum pressure" term in 
the Euler equation [Q. The quantum pressure term is 
important only on length scales of order the coherence 
length, ^ = (87rna)^^/'^, which is typically of order 0.2- 
0.4 /im at the center of the cloud Q|. This is consider- 
ably less than the lengths of the pulses in the experiment, 
which are of order 10 /im, and therefore we shall neglect 
the quantum pressure term in these initial studies. How- 
ever, this term will be important when shocks or other 
sharp structures develop. The equation of motion for a 
is obtained from the continuity equation by integrating 
over the cross section of the cloud, and is 
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The equation of motion for v is found by integrating the 
Euler equation over the cross section of the cloud. 
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where V is the external potential due to the trap and the 
interaction with the laser field. In deriving this result, we 
have assumed that the dependence of on x and y may 
be neglected. We shall assume that the number of parti- 
cles is so large that the Thomas-Fermi approximation is 
valid . For a spatially- uniform, dilute Bose gas at zero 
temperature, the energy per unit volume (and also the 
pressure) are given by E = n^[/o/2. The change in the 
local pressure is thus given by dP — nUodn, where dn is 
the local change in density. The chemical potential, /i, 
neglecting the effects of the trapping potential, is uIIq, 
and in the Thomas- Fermi approximation the sum of this 
and the trap potential V is a constant for a given value of 
z. Thus changes in the chemical potential must likewise 
be functions only of z, and, since dfi = Uodn, it follows 
that changes in the particle density are independent of 
X and y. The change in density is thus given simply in 
terms of the change in the number of particles per unit 
length by dn — da /A, where A is the cross-sectional area 
of the cloud, and therefore the Euler equation becomes 
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Let us begin by considering the linear regime. If we as- 
sume that the variation of the potential in the z-direction 
may be neglected, the problem becomes translationally 
invariant in this direction, and from the linearized ver- 
sions of the continuity equation, Eq. (H), and the Euler 
equation, Eq. (Bj), we find 



d'^a nUo d'^a 
dt^ m dz^ 
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where n = a/Ais the mean density over the cross section 
of the cloud. This result demonstrates that the pulse 
satisfies the wave equation with a velocity, c(z), given by 
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For the case of confinement by a harmonic trap potential 
in the transverse directions, this result has been derived 
independently by Zaremba using different methods. 
In this case n = n{p = 0)/2, where n{p = 0) is the 
density on the axis of the trap. Here p is the radius in 
cylindrical polar coordinates. In the presence of a confin- 
ing potential along the axis of the trap, Eq. (|^) with the 
average density replaced by its local value, n{p = 0, z), 
will be a good approximation, provided the characteris- 
tic length of the disturbance is small compared with the 
axial dimension of the cloud. 

The above derivation underlines the generality of the 
result for the sound speed, but it is instructive to con- 
sider an alternative approach for the specific case of a 



harmonic confining potential in the transverse direction. 
The basic idea is to treat the system as a one-dimensional 
"elastic" medium. The energy of the system consists of 
the kinetic energy due to the bulk motion of the conden- 
sate, an "elastic" energy that takes into account contri- 
butions from particle interactions and the trap potential 
in the transverse direction, and a contribution due to 
the z-component of the trap potential. For simplicity we 
shall assume that the confining potential in the trans- 
verse direction is given by U{p) = where to is 
the frequency of transverse oscillations of a single parti- 
cle in the trap. The kinetic energy, T, per unit length 
is 



T = 



dzi 



dt 
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where C(z) is the displacement of a particle of the fluid 
from its equilibrium position. We next evaluate the elas- 
tic energy, which may be calculated assuming the cross 
section of the cloud to be locally uniform. Per unit 
length, the energy due to the dependence of the trap 
potential on the transverse coordinates is given by 
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Eti-ap = 7:?Ticj / 2t: pdpn{z,p)p , 



(8) 



where the integral is to be taken over the cross section 
of the cloud. The interaction energy per unit volume is 
given by 



^Uq I 2TTpdpn^{z,p). 



(9) 



It is convenient to express quantities in terms of the num- 
ber of particles per unit length. For the harmonic oscil- 
lator trap, the density profile in Thomas-Fermi theory is 
given by [|| 



(10) 



n{p, z) = n(0, z) 1 



R^[z) J ' 



where R{z) is the radius of the cross section of the cloud. 
Thus the number of particles per unit length is given by 



a{z) = in(0,2)7ri?^ 
The density on the axis of the trap is given by 



n(0, z)Uq — -moj^R^, 



(11) 
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and thus one finds that both 7i(0, z) and vary as cr^/^. 
On performing the integrations in Eqs. (^) and one 
arrives at the following expression for the elastic energy 
per unit length: 



E,i = -Uon^i0..z)R' ^ 
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Note that this increases less rapidly with increasing par- 
ticle number per unit length than it would in the case 
of confinement by a rigid pipe with infinitely high walls, 
in which case the potential energy scales as ct^. This 
reflects the fact that for a harmonic confining potential 
an increase in the number of particles per unit length 
results in expansion of the cloud in the radial direction, 
thereby leading to a less rapid increase of the total po- 
tential energy than would be the case for confinement by 
rigid walls. 

Let us now use the result of Eq. ( p^ to calculate the 
velocity of sound in a Bose-condensed cloud when there 
is no potential acting along the axis. The force parallel to 
the axis of the trap is given by F = a^d{Eci/a)/da, and 
the sound speed is given by mc^ = dF/da = ad^Ed/da^, 
results completely analogous to those for the pressure 
and the sound speed in a bulk medium. Inserting the 
expression Eq. (|l^) into this expression, one finds = 
n(p — 0)UQ/(2m). Since for a trap which is harmonic in 
the transverse directions the average density is half the 
density on the axis, this result agrees with Eq. (|). 

We now consider how a pulse will propagate along 
the axis of the trap when allowance is made for the 
varying density. In the Thomas-Fermi approximation, 
the local density on the axis of the trap is given by 
n(0, z) = nniax(l — z'^/Z'^), where nmax is the density 
at the center of the cloud, and Z is the axial dimen- 
sion of the cloud. Thus the local sound speed varies as 
(1 - 2:2/^2^1/2^ ^Yie WKB approximatio n a distur- 
bance initially at z = will at time t have propagated a 
distance z determined by 
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K„axC/o/(2m)]i/2 is 



ox z — Zsin[c(0)t/Z], where c(0) 
the sound velocity, Eq. (0), at z = 0. 

In the experiments relative density disturbances 
were as large as 100%, so it is of interest to investigate 
nonlinear effects. The basic equations are Eqs. (||) and 
(^, which are a closed set of equations for v and cr, 
since from Eqs. ( pi] ) and (12) it follows that the area A 
is given A — 2\Ka\jQl irauP'y^l'^ . We have integrated 
these equations by employing the method of character- 
istics [Q for a number of cases that correspond to the 
experimental conditions. In one set of calculations the 
cloud was initially in equilibrium in a potential due to 
the trap and an extra repulsive potential representing 
the interaction of atoms with the radiation field of the 
laser. The extra potential was taken to have the form 
K-ad = Vo exp(— 2z2/Z2), where / = 12 ^m, as in the 
experiments, and we took the axial extent of the con- 
densate, 2Z, to be 450 /xm. At time f = the repulsive 
potential was turned off. In other runs we started with a 
configuration in equilibrium in the presence of the trap 
alone, and turned the laser field on at i = 0. The strength 
of the potential due to the laser could be varied to pro- 
duce different perturbations in the density. 
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FIG. 1. Density profiles of tlie gas along tlie z-axis for var- 
ious times for cases wlien the perturbing potential is turned 
off, (a), or turned on (b). The profiles are shown for times 0, 
3, 6, 12, 18, 24, 30, 36, 42 and 48 ms. 



In Fig. 1 we show density profiles along the z-axis of 
the cloud for various times for a laser intensity that, in 
equilibrium, would create a relative density perturbation 
of 15% at the center of the cloud. We observe that pos- 
itive density pulses travel faster than negative ones, as 
one would expect from the fact that the sound speed in- 
creases with density. In Fig. 2 we plot the positions of 
extrema of the pulses as functions of time. For reference 
we also show corresponding results for low amplitudes, 
which were obtained by taking the average of the results 
for switching off a potential having the same spatial de- 
pendence, but whose strength was such as to create a rel- 
ative density perturbation of either -1-5%, or —5%. The 
agreement with the result of the linear theory, Eq. (|l4|), 
which is shown as a solid line, is very good, and the small 
difference compared with the numerical results is proba- 
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bly due to the varying cross section of the cloud. Note 
that pulses become distorted as they propagate. This 
is a consequence of nonlinear effects, and the distance a 
pulse can travel before it is distorted significantly may 
be estimated by observing that the sound velocity varies 
as n^/^, where n is the density, and therefore the sound 
velocity change associated with a density change 5n is ap- 
proximately (c/2)5n/n. In a time t a disturbance in the 
region at density n+Sn will move a distance ^ (ct/2)5n/n 
relative to what it would if the density were n, and con- 
sequently distortion of the pulse will be appreciable if 
this distance is greater than or of order the length of the 
pulse, I. In terms of the distance travelled by the pulse, 
z « ci, the condition for there to be significant distortion 
IS z?i2l n/5n^ an estimate which appears to be in accord 
with the numerical calculations. 
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FIG. 2. Distance of the extrema of the density from the 
center of the cloud as a function of time. Squares give results 
for low amplitude pulses, as explained in the text, and the 
solid line is the analytical result, Eq. Diamonds give 

results for the case shown in Fig. la, while crosses correspond 
to that shown in Fig. lb. The central density is taken to be 
lO'^* cm~"^, which corresponds to a time 7r^/[2c(0)] = 96.3 ms 
for a linear pulse to reach the end of the cloud. 



of the particle density over the cross section of the trap, 
and of nonlinearity. Our results indicate that, because 
of their high compressibilities, low density atomic Bose- 
Einstein condensates are useful systems for investigation 
of nonlinear phenomena. 
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FIG. 3. Density profiles of the gas along the z-axis for var- 
ious times for a larger density perturbation. The situations 
correspond to those in Figs, la and lb, apart from the fact 
that the strength of the laser field was chosen so that in equi- 
librium it would create a 30% density perturbation. The pro- 
files are shown for times equal to multiples of 0.85 ms, begin- 
ning at t = 0. 
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To exhibit stronger nonlinear effects, we show in Fig. 3, 
results similar to those in Fig. 1 but for a perturbing po- 
tential that would create a 30% density perturbation in 
equilibrium. The steepening of the leading edge of pos- 
itive density perturbations and of the trailing edge of 
negative ones is apparent. At times later than those for 
which results are shown, shocks develop. To understand 
later stages of the evolution, it is necessary to develop a 
basic understanding of the physics of shocks in superflu- 
ids. In the MIT experiments density perturbations were 
as large as 100%, and under such conditions nonlinear ef- 
fects are expected to be extremely strong. An important 
future experimental challenge is to put in evidence the 
predicted effects on sound propagation of nonuniformity 
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